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Abstract 

Unsustained matter distributions unescapely collapse unless frag- 
mentation and centrifugal or pressure support take place. Starting 
from the above evidence, supermassive compact objects at the cen- 
tre of large-mass galaxies (defined as "holes") are conceived as the 
end-product of the gravitational collapse of local density maxima (de- 
fined as "central collapse") around which positive density perturba- 
tions (overdensities) are located. At the beginning of evolution, as- 
sumed to occur at recombination epoch, local density maxima are 
idealized as homogeneous peaks, while the surrounding envelopes are 
described by a power-law density profile, p{r) oc r^~^, < 6 < 3, 
where b = represents a massless atmosphere and b = 3 a homo- 
geneous layer. The dependence of the density profile on a second 
parameter, chosen to be the ratio between peak and total (truncated) 
mass, K = Mpk/Mtr, is analysed. Overdensity evolution is discussed in 
the context of quintessence cosmological models, which should be use- 
ful in dealing with the virialized phase. Aiming to describe the central 
collapse, further investigation is devoted to a special case where the 
quintessence effect is equivalent to additional curvature {w = —1/3), 
and overdensities exhibit the selected density profile at recombination 
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epoch. A rcdshift-dependent, power-law relation between hole and 
(nonbaryonic) dark halo mass is used to express the dependence of 
the fractional mass, k, on the overdensity mass, M = Mtr, where 
the homogeneous peak and overdensity mass are related to the hole 
and dark halo mass, respectively. Computations are performed for a 
wide range of masses, —1 < log(M/Mio) < 6, and mean overdensity 
heights, 1 < i^i < 4, up to the end of central collapse, and density 
profiles of related configurations arc determined together with addi- 
tional parameters. The central collapse is completed in early times, 
no longer than a few hundredths of Gyr, which implies hole formation 
when proto-haloes, proto-bulges, and proto-disks are still relaxing. No 
appreciable change in evolution (up to the end of central collapse) is 
found with regard to different mean overdensity heights related to 
equal masses. On the other hand, it is recognized that homogeneous 
peaks collapse (in dimensionless coordinates) "faster" with respect to 
surroundings envelopes, in low-mass overdensities than in large-mass 
overdensities. In conclusion, it is inferred that gravitational collapse of 
homogeneous peaks within overdensities may be a viable mechanism 
for hole generation. 

keywords - Dark matter; Dark energy; 95.35.+d; 95.36+x. 

1 Introduction 

The existence of large masses confined in a restricted central region of galax- 
ies, was first suggested in order to explain short-period variability in bright- 
ness exhibited by quasars and active galactic nuclei (e.g., Terrel, 1964; Rees, 
1966), and is widely supported by current high- resolution observations (for a 
review see e.g., Ferrarese and Ford, 2005; Merritt, 2006). There is an amount 
of increasing evidence, that compact objects at the centre of galaxies are su- 
permassive black holes (e.g., Maor, 1998; Miller, 2006), but their origin is 
still debated. 

As a first alternative, supermassive black holes are the end-product of 
compact (high-density) stellar systems which make galactic nuclei, where star 
cncoTintcrs and star collisions arc the dominant physical processes (Spitzer 
and Saslaw, 1966; Spitzer and Stone, 1967; Colgate, 1967; Sanders, 1970). 
Gravitational scattering (via star-star encounters) of energetic stars into elon- 
gated orbits, makes the density of the central region increase which, in turn. 
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implies a higher rate of star evaporation, star colhsions, star disruption, and 
star coalescence. The gas produced during (high-energy) star collisions, star 
disruption, and (high-mass) star death, falls to the central region and con- 
denses into new stars which undergo further collisions. Accordingly, the 
end-product of the evolution of a dense nucleus appears to be the formation 
of a central, massive black hole, cither by runaway star coalescence or grav- 
itational collapse of a massive gas cloud (Begalman and Rees, 1978). The 
possibility of seeds for supermassive black holes which grow via accretion 
of stars and/or gas liberated by star collisions and star disruption, is also 
considered (Duncan and Shapiro, 1983; David et al., 1987a,b; Quinlan and 
Shapiro, 1987, 1989). Seeds for supermassive black holes could also be rem- 
nants of pop. Ill stars (e.g., Madau and Rees, 2001; Johnson and Bromm, 
2007), or the (still speculative) intermediate- mass black holes which form 
in dense star cluster via merger events (e.g., Portegies Zwart et al., 2004; 
Maccarone et al., 2007). For further details, refer to Merritt (2006). 

As a second alternative, supermassive black holes take origin before or 
during galaxy formation. In the spinar model (Morrison, 1969; Cavaliere et 
al., 1969; Cavaliere et al., 1970) the supermassive black hole progenitor is 
a single object (spinar). The spinar may be conceived as a "gravitational 
machine" , where gravitational energy is converted into rotational energy via 
angular momentum conservation during the collapse phase, and rotational 
energy is converted into nonthermal electromagnetic energy (syncroton ra- 
diation) via acceleration of charged particles (mainly electrons) up to the 
relativistic regime, from magnetic torques due to the surface magnetic field. 
According to recent investigations, intermediate-mass black holes acting as 
seeds for supermassive black holes, could be formed from low-angular mo- 
mentum material in primordial disks (Koushiappas et al., 2004) or isothermal 
collapse of atomic hydrogen gas (Mh = lO^-lO^mo) within primordial dark 
matter haloes (Spaan and Silk, 2006). 

The energy released by the formation and growth of supermassive black 
holes must have had a major impact on how gas cooled to form galaxies and 
galaxy clusters (Silk and Rees, 1998). In particular, accretion onto supermas- 
sive black holes provides the energy source for active galactic nuclei which, 
in turn, affects the evolution of galaxies (Silk, 2005). There is increasing evi- 
dence that supermassive black holes play an important role in the formation 
and global evolution of galaxies and of the intergalactic medium (Merloni et 
al., 2005). 
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Seeds for pregalactic black holes could be found much earlier, related to 
topological defects during the inflation epoch, in particular closed domain 
walls (Rubin et al., 2001; Khlopov et al., 2002, 2005; for an extensive review, 
see Khlopov and Rubin, 2004). The resulting mass spectrum extends over a 
wide range, from superheavy to deeply subsolar black holes, where the upper 
limit follows from the condition that pieces of closed walls do not dominate 
within horizon before the whole wall enters it, while the lower limit is given 
by the condition that the black hole gravitational radius exceeds the width 
of contracting domain wall. In addition, low-mass objects are concentrated 
around their most massive counterparts within a cluster (e.g., Khlopov et 
al., 2005). 

Seeds for supermassive black holes, or supermassive black holes, can form 
even in absence of baryons as an inevitable consequence of relativistic core 
collapse following the gravothermal catastrophe, under the restrictive as- 
sumption of self-interacting dark matter haloes (Balberg and Shapiro, 2002). 
The properties exhibited by self-interacting dark matter are in some ways in- 
termediate between their counterparts related to hot and cold dark matter, 
respectively. 

The current attempt starts with the working hypothesis, that supermas- 
sive black holes are the first structures formed within (sufficiently massive) 
density perturbations (hereafter quoted as "overdensities"). Similarly, about 
fifty years ago, galaxies were conceived as emerging from super dense primor- 
dial nuclei (Ambartsumian, 1958, 1965). In fact, unsustained matter distri- 
butions unescapely collapse unless fragmentation and centrifugal or pressure 
support take place. To this respect, the inner central region of collapsing 
proto-galaxies is a very special place, where both tidal effects and acquisi- 
tion of angular momentum are minimized and may safely be thought of as 
negligible. Accordingly, both spherical symmetry and homogeneity are ex- 
pected to be preserved, which implies the formation of a supermassive black 
hole, as in large-mass (m > 20-25mo) supernovae (e.g., Nomoto et al., 2005; 
Ohkube et al., 2006). 

If (nonbaryonic) dark matter is mainly made of fermions (e.g., photinos), 
a fermion ball instead of a black hole could, in principle, take origin (e.g., Vi- 
oUier, 1994; Munyaneza and VioUier, 2002; Munyaneza and Biermann, 2006) 
but no firm conclusion can be reached unless the nature of dark matter is 
known. Concerning ordinary matter, gravitational collapse can be halted 
by Fermi pressure for sufficiently low masses (white dwarfs, neutron stars. 
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quark stars) but the contrary holds for higher masses, where the mean den- 
sity of configurations close to the gravitational i.e. Schwartzschild radius is 
decreased. If a similar trend occurs for dark matter, and supermassive black 
holes are the result of dark matter collapse, a lower mass limit may safely be 
expected. 

A short collapse time of inner, densest regions within ovcrdcnsitics, not 
exceeding a few hundredths of Gyr, would be consistent with the presence 
of a supermassive (M = 3 10^ m©) black hole detected in the quasar SDSSJ 
1148+5251, at a redshift z — 6.41 corresponding to a cosmic age t — 0.840 
Gyr (Willott et al., 2003). 

Supermassive black hole formation via gravitational collapse of the local 
maxima around which overdensities are placed (hereafter quoted as "central 
collapse" ) is a viable scenario, to be tested when the nature and the proper- 
ties of (nonbaryonic) dark matter are known. In this view, the current paper 
investigates the special configuration related to the end of central collapse, 
for different overdcnsity masses (— 1 < log(M/Mio) < 6; Mio = lO^^mo) and 
mean peak heights (1 < I^i < 4), with assigned density profile at the begin- 
ning of evolution. Strictly speaking, a change of the quintessence equation of 
state parameter, and the quintessence degree of clustering, F, should be 
taken into consideration, in dealing with overdensity evolution. The general 
situation shall widely be discussed, to show how the above mentioned pa- 
rameters could affect the virialized configuration. On the other hand, it will 
be shown that central collapse is completed in early times {z > 10) which 
makes the effect of quintessence neghgible concerning both the equation of 
state parameter (e.g., Horellou and Berge, 2005) and the degree of clustering 
(e.g., Nunes and Mota, 2006). Accordingly, detailed calculations shall be 
limited to a special case where considerable simplification may be attained, 
overdensity virialization being outside the aim of the current attempt. 

Some basic considerations are presented in Sect. 2, where the assumed 
initial density profile is also defined. The initial overdensity configuration and 
related evolution up to the end of central collapse, arc described in Sect. 3. 
The results are presented and discussed in Sects. 4 and 5, respectively. The 
conclusion makes the subject of Sect. 6. 
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2 Basic considerations 



2.1 The hole 



The existence of supermassive black holes at the centre of large-mass star 
spheroids, is supported with increasing evidence (e.g., Ferrarese and Ford, 
2005; Merritt, 2006). For sake of brevity, supermassive compact objects 
within a restricted volume (in particular, supermassive black holes) shall be 
quoted as "holes". Accordingly, main components of a typical large-mass 
galaxy are: the (nonbaryonic) dark halo, the stellar halo, the disk, the bulge, 
and the hole. The hole, in turn, is surrounded by an accretion disk, which 
should be quoted as "vortex" or "maelmstrom" , and a dust torus, with the 
emission of a jet along the polar directions during an active phase (e.g., 
Ferrarese and Ford, 2005). 

2.2 Black hole generation via gravitational collapse 

According to the density profile of their progenitor, black hole formation 
can occur bottom-up or top-down. To get more insight, let us refer to a 
spherical-symmetric, truncated power-law mass distribution, as: 



where i?tr is the truncation radius, M^y. the mass within the truncation radius, 
6 = 0, 1, 3 correspond to a Roche sphere (central mass point surrounded by a 
massless atmosphere), a (truncated) isothermal sphere, and a homogeneous 
sphere, respectively. 

Matter distributions expressed by Eq. (1) are simple but nontrivial. For 
instance, models of baryonic collapse within dark matter haloes and for- 
mation of gaseous galactic disks involve use of truncated power-law mass 
distributions (Adams and Bloch, 2006). 

The density profile related to Eq. (1) reads: 




(la) 
(lb) 



0<r<R 



< 6 < 3 




0<r < R, 



(2a) 
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(2b) 



Ptr = 
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where ptr and p^^. are the local density at the truncation radius and the mean 
density within the truncation radius, respectively. The central density cusp 
can be avoided by replacing an arbitrarily small inner core at r = Tpk, with 
a homogeneous core where p{r) = p(rpk), < r < Tpk- 

A spherical-symmetric matter distribution collapses into a black hole 
when the gravitational i.e. Schwartzschild radius (e.g., Landau and Lifchitz, 
1966, Chap. XI, §97) is attained, as: 



where G is the constant of gravitation and c the light velocity in vacuum, to 
be intended as the light velocity in absence of baryonic matter. For further 
details, refer to Appendix A. 

For the density profiles under discussion, the gravitational radius related 
to a sphere of radius, r, enclosing a mass, M(r), is obtained by the combina- 
tion of Eqs. (1) and (3), the latter generalized to a mass, M(r). The result 
is: 



where Rg^ is the gravitational radius related to the truncation mass, according 
to Eq. (3). 

In the special case of the isothermal sphere, 6 = 1, all the spheres collapse 
into a black hole at the same time. Milder density profiles, 1 < 6 < 3, imply 
Tgr/r < Rgr/Rtr i-G. top-down black hole formation: the gravitational radius 
is first attained at the truncation radius. Harder density profiles, < 6 < 1, 
imply Tgr/r > i?gj./-Rtr i-e. bottom-up black hole formation: the gravitational 
radius is first attained at the centre. 

If holes are the end-products of overdensity central collapse, the formation 
mechanism must necessarily be bottom-up with regard to the whole mass, but 
all the possibilities remain open with regard to the local maximum, conceived 
as the proto-hole. 

2.3 Initial density profiles 

Power-law density profiles provide a poor fit to dark matter haloes unless 
sufficiently thin shells are considered. A better description relies on three- 
exponent laws with a central cusp (e.g., Hernquist, 1990; Navarro et al.. 




(4) 
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1995, 1996; Zhao, 1996; Moore et al., 1998, 1999; Fukushige and Makino, 
2001, 2003; Diemand et al., 2004; Reed et al., 2005; for the description of a 
fitting procedure see e.g., Caimmi and Marmo, 2003, 2004; Caimmi et al., 
2005; Caimmi, 2006) or two-parameter laws (e.g., Sersic, 1968; Navarro et 
al., 2004; Merritt et al., 2005). 

On the other hand, it is not the case at the beginning of the evolution 
which, for high-energy universes, may safely be related to recombination 
epoch (e.g., Caimmi, 1989), where the contribution of dark energy is still 
negligible. Aiming to describe overdensity evolution up to the end of the 
central collapse, a simple but nontrivial density profile is assumed. 

More specifically, overdensities at recombination epoch {z ^ 1100) are 
represented as homogeneous peaks surrounded by envelopes where the den- 
sity profile is represented by a power- law, according to Eqs. (2). The exphcit 
expression reads: 



p(r) 



6-3 



< r < Tpk ; 
J^pk < r < Rtr ; 



(5) 



where and are the homogeneous peak density and radius, respectively. 

Owing to spherical symmetry, the mass enclosed within the truncation 
surface is: 

Mtr = / p(r)47rr^ dr ; (6) 
Jo 

which may be integrated using Eq. (5). The result is: 



Mtr = Mpk <! 1 + 

4:71 

Mpk = yPpkrpk ; 




- 1 



(7a) 
(7b) 



where Mpk is the homogeneous peak mass. 

Using Eq. (7a), a dimensionless radius, S, may be defined in terms of a 
dimensionless mass, as: 



b l-K 

3 K, 



1/6 



M 



K 



pk 



M 



(8a) 
(8b) 
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and the mean density within the truncation radius, p^j. = (3/47r)(Mtr/i?tr)) 
owing to Eqs. (7) and (8), takes the expression: 



in terms of the dimensionless radius, S, and: 



1 + 



bl-K 

3 K 



-3/b 



(9) 



(10) 



in terms of the dimensionless mass, k. 

The combination of Eqs. (5) and (8) yields the local density at the trun- 
cation radius, as: 

Ptr 



Ppk 



bl-K 

1 + 7T 

6 K 



{b-3)/b 



11) 



and the combination of Eqs. (10) and (11) yields the local to global density 
ratio at the truncation radius, as: 



Ptr 



K 



3 K 



(12) 



which depend on two parameters, b and k. 

With regard to a generic dimensionless radius, ^, .^pk < C < ^, Eqs. (7), 
(9), and (5) may be generalized as: 



M(r) = Mpk 

-Pir) - f 
p{r) = PpkC 



6-3 



and the combination of Eqs. (14) and (15) yields: 

n -1 



p(r) 
p(r) 



3 6-3^ 
6 + ^^" 



(13) 
(14) 
(15) 

(16) 



which is the generalization of Eq. (12) to a generic dimensionless radius, ^, 
U < ^ < 

For a complete definition of density profiles within overdensities at re- 
combination epoch, the boundary conditions derived from the cosmological 
model are needed. 
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3 Overdensity evolution up to the end of cen- 
tral collapse 

Let overdensities be made of a homogeneous peak surrounded by a heteroge- 
neous envelope, at recombination epoch. The vahdity of the above mentioned 
assumption imphes hole formation unless the collapse is halted by Fermi pres- 
sure, as in neutron star progenitors. In absence of any information about the 
(nonbaryonic) dark matter, hole formation via gravitational collapse cannot 
be escluded and must be taken into consideration. To this aim, overden- 
sity evolution up to the end of central collapse shall be investigated in what 
follows. 

Given a spherical-symmetric overdensity around a local maximum at the 
beginning of evolution, the central collapse is intended as the collapse of a 
dusty homogeneous sphere with same density and expansion rate with respect 
to the local maximum. 

3.1 Overdensities at recombination epoch 

The effect of dark energy may safely be neglected in early times {z > 10) 
and matter may be described, to a good extent after recombination, in the 
context of Fricdmann universes with dust only. With regard to overdensities 
bound around a local density maximum, or peak (e.g.. Heavens and Peacock, 
1998), related mean and local values are defined as: 

; 5^-^-1 ; (17) 

Phm Phm 

where is the mean density of the matter Hubble flow. Keeping in mind 
the definition of mean density, the combination of Eqs. (17) yields (e.g., Ry- 
den and Gunn, 1987): 

where S is the overdensity volume. 

For high-energy universes, as in current QCDM cosmologies where flm ~ 
0.3, Qq ~ 0.7, Qm + ^5 = Ij the initial configuration may safely be related 
to the recombination epoch (see e.g., Caimmi, 1989 for the special case of 
CDM cosmologies). According to Eq. (17), the overdensity mass reads: 

M ={l + 6i)phm{ahi)Si = {l + 6i)Mhm ; (19) 
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where Mhm is the mass of matter Hubble flow enclosed within the volume, 
Si, and the index, i, denotes recombination. 

The mean overdensity value can be expressed in terms of the peak over- 
density value, 5pk, as: 



5 = F5, 



pk ) 



1 f f f Sjr) 
S J J Js (5pk 



(20) 



where the shape factor, F, lies in the range, Sp]^/ S < F < 1, between the 
extreme situations of a peak surrounded by a masslcss atmosphere (lower 
limit) and by a layer of equal density (upper limit); and S'pk represents the 
volume filled by the peak, assumed to be homogeneous. 

Let the mean and the local overdensity height be defined as the ratios: 

^ ^ SM-<f>li' ; (21) 



6 



M 



'M 



— 2 1/2 

where < 6 >^,i is the rms overdensity value related to an assigned volume, 
S, at a given epoch (recombination in the case under discussion), and the 
index, M, means that a volume, S, within the Hubble fiow, encloses a mass, 
M = Mhm = PhmS- The combination of Eqs. (20) and (21) yields: 

i^ = Fi^p^ ; i/pk = fc . (22) 

If, in addition, the density profile is expressed by Eq. (5), the combination 
of Eqs. (14), (15), (16), and (17) yields: 



1 + 5(r) 



1 + 6, 



pk 



l + 5(r)-(l + 5pk)e'-' 
'3 



1 + 6{r) = [1 + 6{r) 



Ppk — Phm 

Phm 



6-3, 
6 + ^^" 



(23) 
(24) 

(25) 
(26) 



and the particularization to the truncation radius, ^ = S, by use of Eqs. (7a) 
and (8) produces: 

-3/b 



1 + 5tr = 



1 + S- 



pk 



1 + 



b 1 



K 



(27) 
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1 + 5tr = (1 + V) 



l + 5tr = {l + 6tr)K 



bl-K 

6 K 

bl-K 



1-3/6 



1 + 



3 K 



lim 

6^0 



< 



^ bl-K 

3 K, 

bl-K 

1 + 



-1 1-3/6 



= exp 



1- K 



K 



3 K 



1-3/6 



<i ; 



< 6 < 3 



(28) 
(29) 
(30) 
(31) 



accordingly, overdensities at recombination are defined by the following pa- 
rameters: total mass, M = Mti-; mean value, 6 = Str', mean height, V = V^r'-, 
fractional mass, k] power-law density profile exponent, 6; which allow the 
calculation of the peak overdensity, 5pk, and the local overdensity, 5(r). At 
this stage, a cosmological model needs to be specified. 



3.2 Input parameters and cosmological model 

Overdensity masses and mean heights shall be considered in the range, — 1 < 
log(M/Mio) < 6 (Mio = lO^^mg) and 1 < V < 4, respectively. Present-day 
rms overdensity values in the linear approximation, shall be assumed as in 
CDM cosmologies (e.g., Gunn, 1987; Ryden and Gunn, 1987). Related values 
at recombination may be determined, provided the temporal behaviour of the 
rms overdensity value for a selected mass, is known. 

For QCDM cosmologies, it has still to be found an analytical solution 
which represents the growing mode of overdensities. In the special case of 
flat universes, Qm + = 1, and time-independent quintessence equation of 
state parameter, w = const, the solution may be written in terms of the 
hypergeometric function, 2F1 (Silveria and Wega, 1994). 

Extending the results related to the special case where the dark energy 
mimics a cosmological constant, w = —1 (Wang and Steinhardt, 1998; Carroll 
et al., 1992) to the range, w > —1 (Basilakos, 2003), and improving the 
results related to phantom energy, w < —1, yield the following approximation 
(Percival, 2005): 

Sujah) _ 5 Qrn{ah)ah . ^ 

^m(I) ~ 2 [f)„(aO]" - ^,{ah) + [1 + 0.5f7^(aO] [1 + X^M] ' ^ 
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Q f{o-h) 

= / [1 + w{x)]dHx) ; (32d) 

3(1 - w) f 1 2-3w , ^ , ,,1 , , 

_ 1 75.. + 76 . 

n^{ah) + ^q{ah) ^ 1 ; n„ = n^(l) ; ^5 = 0^(1) ; (32g) 
w{ah) — const = ty ; (32h) 

where ah is the cosmological scale factor, normahzed to = 1 at the present 
time, and are density parameters related to matter (including radia- 
tion in the case under discussion, z ~ 1100) and quintessence, respectively. 
In the limit of constant quintessence equation of state parameter, Eq. (32d) 
reduces to (e.g., Horellou and Berge, 2005): 

/(a^) = -3(l + ^) ; (33) 

according to Eq. (32h). 

At recombination epoch, z — 1100 to a good extent, and the cosmological 
scale factor reads: 

where the index, i, denotes the recombination. 

The input parameters of the cosmological model are assumed to be the 
following (e.g., Manera and Mota, 2006): 

fl^ = 0.3 ; flq = 0.7 ; fl^ = 0.047 ; h = 0.65 ; cxg = 0.9 ; 

(35) 

where Qi, is the baryon density parameter, h — if/(100kms ^Mpc ^) is the 
dimensionless Hubble parameter at 2; = 0, and as = Sm{8) is related to the 
mass, M(8), within a sphere of radius, R{8) = 8/i~^ Mpc, filled by Hubble 
flow at present. 
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A different dependence, ag = f{w), is deduced according if X-ray cluster 
data (Wang and Steinliardt, 1998) or cosmic microwave background data 
(Dor an et al., 2001) are used. The cases, ag — 0.8 and erg = 0.7 should also 
be considered (e.g., Horellou and Berge, 2005), but a complete investigation 
on this point is outside the aim of the current paper. 

The system of measure, [kpc Mio Gyr], Mio = 10^'^mQ, shall be adopted 
in performing computations. The conversion formulae from and towards the 
standard astrophysical system of measure, [cm g s], are reported in Appendix 
B. 

Using Eqs. (35), the current value of critical and mean matter density are: 

Pcrit = ^ = 1.172 6 10-8Miokpc-3 ; (36a) 
P/^m = ^^mPcrit = 3.517810-9Miokpc-' ; (36b) 

and the amount of matter Hubble flow enclosed within a sphere of radius, 
R{8) = 8/i-^Mpc, is: 

47r / 8 \ ^ 

M{8)^—Phm[j^) = 94402Mio ; (37) 

which is consistent with a value M(8) = 932513Mio deduced by use of a 
CDM cosmological model (Gunn, 1987). The rms overdensity spectrum at 
recombination epoch, 5m (oh?), shall be deduced from the above mentioned 
GDM model (e.g., Gaimmi et al., 1990). Related present-day values in the 
linear theory, 5m(1), are deduced from the standard relation (e.g., Zeldovich 
and Novikov, 1982): 

Phmiah) ^ Phmil)il + zf ; (38) 

where z is the redshift. 

The radius of a sphere filled by matter Hubble flow of mass, M — lO^Mio, 
is defined by the relation: 

4:71 

M=10^Mio = p/.^(a,,)— i?3(a0 ; (39) 

and the combination of Eqs. (36) and (39) yields: 

Rq = R(l) = 0.189 327 10(^+^°)/^kpc ; (40) 
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with regard to the current epoch, where —1 < £ < 6 in the range of interest. 
A hnear dependence, R oc ah, imphes the following: 

i^M-jf^ ; (41) 

at any selected time; in particular, Ri = i?o/1101 at recombination epoch. 

For nonrotating, spherical-symmetric overdensities with null peculiar ve- 
locity field, the turnaround radius, rmax, is defined by the relation (e.g., 
Peebles, 1980, Chap. II, § 19A; Caimmi, 1989): 

^ = A, ; (42a) 



A, = l-((]^)ri(l + 5,)-i>0 ; (42b) 
which holds for CDM universes. The combination of Eqs. (21) and (42) yields: 

= i?ipMi)~' ; |l-(^]„^)-^|<^^<l ; (43) 



to the first order in Si. 

The rms overdensity value in the linear theory, Sui, and the radius of the 
related sphere, R, are listed in Tab. 1 with regard to current and recombina- 
tion epoch, using a CDM model (Gunn, 1987), as a function of the mass. The 
product, r'^rmax, deduced from Eq. (43), is also listed therein. Overdensity 
values at recombination epoch, Sui, deduced from their present-day coun- 
terparts, 6mo, by use of Eqs. (32), related to QCDM models with constant 
quintessence equation of state parameter, w, yield Zi = 1101.5263 instead of 
the assumed value, Zi = 1100, making a close agreement, as expected. 

Turning to the general case of QCDM cosmologies, the dynamical expan- 
sion of the universe is described by the Priedmann equations (e.g., Percival, 
2005): 



Qfc = 1 - fim - fig ; ^u{ah) = ; u^m,q ; (44b) 

Pcrit[0'h) 

PcHt(a.) = ^^^ ; H{ah)^^ ; (44c) 
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log(M/Mio) 




i?o/kpc 






i?i/kpc 




-1 


14 


1.89 E+2 


1.27 E- 


-2 


1.72 E-1 


8.13 E+0 





11 


4.08 E+2 


1.00 E- 


-2 


3.70 E-1 


2.23 E+1 


1 


8.0 


8.79 E+2 


7.26 E- 


-3 


7.98 E-1 


6.60 E+1 


2 


5.3 


1.89 E+3 


4.81 E- 


-3 


1.72 E+0 


2.15 E+2 


3 


3.4 


4.08 E+3 


3.08 E- 


-3 


3.70 E+0 


7.21 E+2 


4 


1.9 


8.79 E+3 


1.72 E- 


-3 


7.98 E+0 


2.78 E+3 


5 


0.87 


1.89 E+4 


7.89 E- 


-4 


1.72 E+1 


1.31 E+4 


6 


0.32 


4.08 E+4 


2.90 E- 


-4 


3.70 E+1 


7.66 E+4 



Table 1: The rms overdcnsity value in the linear theory, 6u, and the radius 
of the related sphere, R, with regard to current and recombination epoch, 
using a CDM model (Gunn, 1987), as a function of the mass. The product, 
T^i^max, deduced from Eq. (43), is also listed. The current and recombination 
epoch are denoted by the indices, and i, respectively. Overdensity values 
at recombination epoch, Sui, deduced from their present-day counterparts, 
^MOi by use of Eqs. (32), related to QCDM models with constant quintessence 
equation of state parameter, w, yield Zi = 1101.526 3 instead of the assumed 
value, Zi — 1100, making a close agreement, as expected. 
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where fi^ is the curvature density parameter, which is null for flat universes, 
and Eqs. (32) hold. The function, /(a/j), is calculated by solving the conser- 
vation of energy equation for the dark energy (Caldwell et al., 1998a): 

= -3pA ; (45) 

which yields: 

p,(a,) = pX^"") ; (46) 
on the other hand, the conservation of matter mass produces: 

Pm{ah) = Pm^l^ ; (47) 

in the case under consideration of flat universes, the combination of Eqs. (44), 
(46), and (47) yields: 

© ' ^ ^'^"'^^ ^ ^ K"'^' + ^^"^'"^] ' ^^^^ 

where the function, /(o/j), is defined by Eq. (32d), or Eq. (33) in the limit of 
constant quintessence equation of state parameter, w. 

Though the action of the dark energy may be explained in different 
ways (e.g., Caldwell et al., 1998a), the current interpretation is based on 
a dynamical evolving scalar field slowly roUing down its potential, V(0). If 
quintessence is minimally coupled to gravity i.e. any couplings to other fields 
are supposed to be negligibly small, the equation of motion for the scalar field 
within the Hubble flow is (e.g., Manera and Mota, 2006): 

dV 

^=-^H4>-— ; (49) 

where the variable, 0, may be conceived as a curvilinear coordinate (with 
the dimension of a length) which describes the evolution of the scalar field. 
A dynamical analogon of Eq. (49) is e.g., a solid shpere in free fall within a 
homogeneous fluid under the action of gravitation, where the flrst and the 
second term on the right-hand side of Eq. (49) are related to viscous and 
gravitational force, respectively. 
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The pressure and the density of the scalar field are (e.g., Nunes and Mota, 
2006): 

= ^PK0' - ppV(0) ; (50) 

P. = ^f^^ + fv(0) ; (51) 

where px and pp are mass densities related to kinetic and potential energy, 
respectively, which implies pk0^/2 and ppV{(f>) are kinetic and potential en- 
ergy densities, respectively, of the scalar field^. It is apparent that a slow 
rolling of the scalar field, pk0^/2 <C ppV{4>), implies a negative pressure via 
Eq. (50). 

The quintessence equation of state parameter, via Eqs. (50) and (51), 
reads: 

w^^^ PK0V2 - ppV(0) _ 
c^p, pK0V2 + PpV(0) ' 

where, in general, w changes in time. In the special case of a static scalar 
field, = 0, w = — 1, and the quintessence mimics the effect of a cosmolog- 
ical constant. In the special case of a steady rolhng, =const, and time- 
independent potential, V(0) = (1/s)(pk/pp)(</'/2), s real number, Eq. (52) 
reduces to: w = {s - l)/(s + I), where s = 0, 1/6, 1/5, 1/4, 1/3, 1/2, 1, 
yields w = —1, —5/7, —2/3, —3/5, —1/2, —1/3, 0, respectively, w is also 
time-independent, and the quintessence mimics the effect of an ecsessence 
(e.g., Iliev and Shapiro, 2001; Horellou and Berge, 2005), which is a perfect 
fluid (e.g., Nunes and Mota, 2006). 

The general case, w — w{t), can result from a changing ratio of quintessence 
kinetic to potential energy which, in turn, is owing to the evolution of the 
related scalar field potential (e.g., Caldwell et al., 1998a). Though time- 
varying equations of state are closer to the real situation (e.g., Wetterich, 
1995; Amendola, 2000; Battye and Weller, 2003; Mota and van de Bruck, 
2004; Percival, 2005; Manera and Mota, 2006; Nunes and Mota, 2006; Maio 

^The mass densities, pK and pp, are usually omitted in literature under the pretext 
that the chosen units imply unitary values of the light velocity in vacuum, c = luc, and 
the Planck constant, h = luh, where Uc and uu are the unit velocity and the unit action, 
respectively, related to the chosen system of measure. In author's opinion, this practice 
is extremely dangerous, as the dimensions of the terms appearing in any equation cannot 
be tested. 
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et al., 2006; Basilakos and Voglis, 2007), the special case of constant w makes 
considerable simplification (e.g., Caldwell et al., 1998a; Wang and Steinhardt, 
1998; Weinberg and Kamionkowski, 2003; Horellou and Berge, 2005; Maor 
and Lahav, 2005). Qualitatively, many of tfie effects of general w{t) models 
can be predicted by interpolating between models with constant w. If the 
quintessence equation of state varies only slowly with time, the observational 
predictions are well approximated by treating w{i) = const = w (Wang et 
al., 2000; Percival, 2005). More realistic models where w — w{t) such as 
2EXP (e.g., Nunes and Mota, 2006) and SUGRA (e.g., Maio et al, 2006) 
show that w = const to a good extent for z ~ 10 or i ~ 0.5 Gyr. 

In dealing with overdensity central collapse, time scales no longer than 
a few hudredths of Gyr are involved, or redshift substantially larger than 
about z — 10. Accordingly, w — const may safely be assumed in the current 
paper. 

The combination of Eqs. (33), (46), (47) and (48) yields at recombination: 



where Rh is the radius of a sphere filled by Hubble flow of mass, M^h + 
Mghio-h), and Mg is formally defined as a quintessence "mass" (Caimmi, 
2007). Let a few limiting but relevant situations be illustrated with more 
detail. 

In the special case, w — Q, f{ah) = —3, according to Eq. (33), and the 
effect of quintessence is equivalent to the presence of additional matter, as 
shown by Eq. (44a). Nothing changes with respect to a CDM universe where 
(f2m)cDM = + (e.g., Caldwell et al., 1998a). The quintessence mass is 
time-independent via Eq. (53b). 

In the special case, w — —1/3, f{ah) = —2, according to Eq. (33), and the 
effect of quintessence is equivalent to the presence of additional curvature, as 
shown by Eq. (44a). Nothing changes with respect to a CDM universe where 
(^fc)cDM = ^k + ^q (e.g., Horellou and Berge, 2005). The quintessence mass 
scales as ah via Eq. (53b). 

In the special case, w = —1, f{ah) = 0, according to Eq. (33), and the 
effect of quintessence is equivalent to the presence of a cosmological constant. 





(53b) 



(53a) 
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as shown by Eq. (44a). Nothing changes with respect to a CDM universe in 
presence of a cosmological constant where A/ (3if ^) = Qq (e.g., Horellou and 
Berge, 2005). The quintessence mass scales as via Eq. (53b). 

3.3 Overdensity expansion and central collapse: gen- 
eral ideas 

The dynamics of nonhnear structure formation in the universe may show a 
distinct signature associated to the nature of the dark energy and a partic- 
ular model (Mota and van de Bruck, 2004). In this view, the behaviour of 
quintessence during the nonlinear regime of structure formation, is conceived 
as lying between two limiting cases, namely (i) full clustering i.e. the scalar 
field responds to the infall in the same way as the matter, and (ii) unclus- 
tering i.e. the scalar field remains homogeneous and the sole effect is a tidal 
potential acting on the matter overdensity. 

The general case of partial clustering can be taken into consideration, 
at the expense of a much more complicated continuity equation for the 
quintessence overdensity (e.g., Mota and van de Bruck, 2004; Maor and La- 
hav, 2005; Nunes and Mota, 2006). Even though it has been shown that 
quintessence cannot be perfectly smooth (Caldwell et al., 1998a,b), cluster- 
ing is usually assumed to be negligible on scales less than about 100 Mpc 
(e.g., Wang and Steinhardt, 1998; Weinberg and Kamionkowski, 2003; Bat- 
tye and Weller, 2003; Horellou and Berge, 2005). It is therefore common 
practice to keep the quintessence homogeneous during the evolution of over- 
densities. The effects of relaxing this assumption were explored in recent 
attempts (e.g., Mota and van de Bruck, 2004; Percival, 2005; Manera and 
Mota, 2006; Nunes and Mota, 2006). 

Given a spherical-symmetric overdensity and an infinitely thin spheri- 
cal shell at a distance, r, from the centre, in the limit of fully clustered 
quintessence, the equation of motion is still expressed by Eq. (48) provided 
the Hubble parameter, H, is replaced hj Hr = r /r (e.g., Nunes and Mota, 
2006). In general, the evolution of the density of the scalar field (continuity 
equation) within the shell is (e.g., Mota and van de Bruck, 2004; Maor and 
Lahav, 2005; Nunes and Mota, 2006): 

Pg + 3-{l + w)pg^r ; (54) 
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where T describes the quintessence density change due to the cosmic expan- 
sion. In the limit of fully clustered quintessence, F = as for matter. In 
the limit of unclustered quintessence, F = —3{ah/ah — r/r)(l + w)pg, and 
Eq. (54) reduces to its counterpart related to the quintessence Hubble flow. 

The shell expansion is related to cosmic expansion and overdensity change, 
as (e.g., Nunes and Mota, 2006): 



where 6{r) is the shell overdensity, defined by Eq. (17). 

The real situation may safely be expected to lie between the above men- 
tioned limiting cases, where the quintessence is clustering together with the 
matter and remains homogeneous, respectively (e.g., Maor and Lahav, 2005; 
Caimmi, 2007). In general, ovcrdensities reach turnaround and collapse ear- 
lier in QCDM models with larger quintessence equation of state parameter, 
w, and vice versa (e.g., Weinberg and Kamionkowski, 2003; Horellou and 
Berge, 2005). On the other hand, the difference is small in early times, 
where the effect of the dark energy is still negligible (e.g., Weinberg and 
Kamionkowski, 2003; Horellou and Berge, 2005; Nuncs and Mota, 2006). 

For this reason, the quintessence equation of state parameter, w, shall be 
kept constant during overdensity evolution up to the end of central collapse, 
which is expected to occur in early times. The above condition, in turn, im- 
plies that different evolutions related to different degrees of quintessence are 
essentially indistinguishable except for the later stages when quintessence 
starts to dominate at low redshifts. In general, the evolution depends on 
the quintessence equation of state, the potential of the scalar field, and the 
quintessence contribution to the total energy budget of the universe at high 
redshifts. For further details, refer to e.g., Nunes and Mota (2006). Ac- 
cordingly, the quintessence within ovcrdensities shall be assumed as fully 
clustered to simplify the calculations. 

Under the above mentioned restrictions, the density ratio of clustered to 
unclustered quintessence may be expressed in terms of the related overdensity 
value, as (e.g., Nunes and Mota, 2006): 



r dh I S{r 




(55) 




w — const 



(56) 
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where p^/, is the density within the quintessence Hubble flow and the index, 
i, denotes the beginning of evolution. 



The integration of Eq. (54) where w — const, F = 0, yields (e.g., Caimmi, 
2007): 



with regard to the turnaround configuration, r = rmax- Accordingly, the 
quintessence density is related to the radius by a power law with exponent, 
b — 3 — —3w — 3. The special case, w — —1/3, b — 1, corresponds to an 
isothermal sphere; w — —1, 6 = 3, to a homogeneous sphere; w — 0, b — 0, 
to a Roche sphere. 

The quintessence mass, Mq{r), enclosed by a homogeneous spherical over- 
density with equal density and radius as in the infinitely thin sperical shell 
under consideration, scales as: 



which makes the mass increase as the radius increases, and vice versa. 

The validity of Newton's theorem and McLaurin's theorem being inde- 
pendent of the value of the gravitation constant, G, homogeneous overden- 
sities where the quintessence is fully clustered evolve in the same way as 
within static quintessence density profiles, defined by Eqs. (57) or (58), even 
if the value of the interaction strength is different. In this view, Mg{r) is the 
quintessence mass related to the static profile, enclosed within the volume of 
the homogeneous over density of radius, r. 

Owing to the above mentioned theorems, (i) given a quintessence spherical 
corona, related to the static profile, the gravitational action on a selected 
point enclosed by the inner surface of the corona, is null and (ii) given a 
quintessence sphere related to the static profile, enclosing a mass, Mg(r'), 
the gravitational action on a selected point outside the surface is equivalent 
to its counterpart exerted by a central point of equal mass, Mq(r'). 

Using the above results, the evolution of overdensity infinitely thin, sh- 
perical shells may be determined in the limit of fully clustered quintessence. 
To this aim, the gravitational potential induced by the quintessence is con- 
ceived as arising from a distribution where any two "particles" , idealized as 




w — const 



r = o 



(57) 




(58) 
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"mass points", rriqi and rrigj, interact with strength, (1 + 3w)G, according 
to a Newton-hke law, Fij = (1 + 3w)G'mqimqj/rij. Then the results related 
to two-component matter distributions (e.g.. Limber, 1959; Brosche et al., 
1983; Caimmi et al., 1984; Caimmi and Secco, 1992) may be generalized to 
the case, where a subsystem is made of quintessence. For further details, 
refer to Caimmi (2007). 

In the special case under consideration {w = const, F = 0) the expansion 
of an infinitely thin spherical shell enclosing a matter mass, Mjn{r), and a 
quintessence mass, Mg(r), within a radius, r, is expressed as (e.g., Mota and 
van de Bruck, 2004): 

_^gjU0^(l + 3^)GM.W^^ , (69) 

where mass conservation holds for matter, M^^^{r) = Mm{ri), and quintessence 
mass changes in time according to Eq. (58). The special cases, w — and 
w — —1/3, relate to flat (JIcdm = + = 1) and open (JIcdm = — 
1 — Qq) CDM universes. 

The quintessence mass, Mq{r), is assumed to change in time according 
to Eq. (58), which is related to the quintessence mass enclosed by a homoge- 
neous spherical overdensity with equal density and radius as in the infinitely 
thin shell under consideration. The above approximation overstimates the 
collapse rate, where the larger effect corresponds to the lower quintessence 
equation of state parameter, and vice versa. In any case, the effect of the dark 
energy is negligible in early times (e.g., Weinberg and Kamionkowski, 2003; 
Horellou and Berge, 2005; Nunes and Mota, 2006), when central collapse is 
completed. 

Following a similar procedure as in presence of sole matter (e.g., Caimmi, 
1989), an integration of Eq. (59) yields: 

r^Tin \ 

[ n r 

2(1 + 3«;)GM,(r^a.)(rmax)'" , 2(1 + 3w)GMq{r^,^,)(r^^f^ \ 



where the plus is related to expansion, and the minus to contraction. 

To a good extent, at recombination epoch overdensities still expand to- 
gether with the universe, which implies the boundary conditions: 

U = Thi ; u = fhi ; Thi = auTho ; (61) 
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where ry, is the radius of a sphere with zero overdensity. 

The particularization of Eq. (48) to au = dhi, after combination with 
Eqs.(19), (33), (46), (47), (58), and (61) yields: 



' max y V max J 



3w 



1/2 



(62) 



provided isotropic radial motions arc dominant. 

At turnaround, r = by definition, which implies by use of Eq. (60) 
particularized to r = rjnax' 



.2 2GM^{n) 2(l + 3«;)G'M,(w)( 

' max J 



2GM^{n) 2(1 + 3w)G'M,(r^ax) 

'"max '"max 

and the substitution of r? into the general expression, Eq. (60), yields: 
^2GMr^{n) 2GM^{n) 



(63) 



2(1 + 3^/;)G'M,(r^ax)(rmax)'" 2(1 + 3«;)G'Mg(r„ 



j.l + 'iw 



1/2 



; (64) 



where the turnaround radius, rmax, has replaced the recombination radius. 
Using the dimensionless variables: 



t 



a — 



T — 



8G'M^(r,) 



tmax 
.-1/2 
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-1/2 



(65a) 
(65b) 



together with the fractional mass: 



Mg(rmax) 



m = m(^rmaxj 

1 Mm(rj) m Mq{r 

max > 



Mmiri) ' 



1 + m Mtot {r max 

) 1 + m Mtot 



(66a) 

■; M^^t^Mm + Mg; (66b) 
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allow to cast Eq. (64) under the equivalent form: 



da TT 
d7 " ^2 



i-l)+(l + 3»)m(-^-l) 



1/2 

; (67) 



where the special cases, w = and w = —1/3. reproduce the results related 
to flat (i^cDM = + = 1) and open (r2cDM = — ^ — ^q) CDM 
universes (e.g., Ryden and Gunn, 1987; Caimmi, 1989). 

It is worth noticing the scaling time, t]^^, defined by Eq. (65b), coincides 
with the actual matter overdensity free-fall time only for CDM universes 
i.e. Mtot(rmax) = -Hn(rmax), Mg(rmax) = 0. In the case under discussion, 
overdensity collapse is equivalent to pure matter collapse within a static 
quintessence density profile, Pg(r)/pg(rmax) = {r / rraax)~^^^~^^\ according to 
Eq. (57). The related free-fall time, to be calculated numerically, is the over- 
density free-fall time. Then the dimensionless time, r, is different from unity 
at turnaround, unless CDM universes are considered (e.g., Ryden and Gunn, 
1987; Caimmi, 1989). In any case, CDM models make a useful reference case, 
and for this reason the main results shall be recalled in the following. 

In the special case, w — 0, the effect of quintessence is equivalent to the 
effect of additional matter, JIcdm = + = 1, and after a re-definition 
of the dimensionless time, Eq. (67) reduces to: 

da 71 \J OL — oP- . 

r = Vl + m^ ; (68b) 

where the scaling time, t\a^l Vl + ^5 represents the overdensity free-fall time 
with respect to the turnaround configuration. An integration yields (e.g., 
Caimmi, 1989): 



2 

T - Tj = ^- 
TT 



-\J a — o? -|- ^Jai — of — arcsin Vl — cc -|- arcsin Vl — Oii ; 

(69) 

in the idealized situation where radial motion is preserved, a — 1 and r = 
2i — 1 at i-th turnaround; a = and r = 2i at i-th point-like configuration; 
and Eq. (69) can be splitted as: 



TT 



T = (2i — 1) Va — a^ -l- arcsin — a ; (70a) 
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T — {2i — 1) -\ — — + arcsin ^/T 



a 



TT 



(70b) 



which describe i-th expansion and i-th contraction, respectively. 

For sake of completeness, it is worth recalling the parametric equations 
(e.g., Ryden and Gunn, 1987; Caimmi, 1989): 



^-sin^ 



1 — cos 9 



T — 

dr _ 2 

d^ ~ TT 



a — 



TT 



a 



da 
"d^ 



(71a) 
(71b) 



where Eq. (68) has been used for the expression of da/ d^. 

In the special case, w = —1/3, the effect of quintessence is equivalent to 
the effect of additional curvature, ficDM = = 1 — ^q, and Eq. (67) reduces 
to: 



da 

d7 



TT 



a — a 



(72) 

2a' ^ ' 

which coincides with Eq. (68a), and the validity of Eqs. (69)-(71) is main- 
tained. 

In the special case, w = —1, the effect of quintessence is equivalent to the 
effect of a cosmological constant. A/ (SiJ^) = fig, and Eq. (67) reduces to: 



da TT 
~ ^2 LVa 



--l)-2m(a^-l) 



1/2 



(73) 



which has no analytical solution and must be solved via numerical integra- 
tion. 



3.4 Overdensity expansion and central collapse: a spe- 
cial case 

To the aim of the current attempt, only the early phase of overdensity evolu- 
tion is considered, from recombination to the end of central collapse, which 
coincides with the onset of shell crossing in the limit of classical mechanics, 
where black hole formation does not occur and radial motions are reversed 
off from the centre. The central collapse may safely be expected to end in 
a considerably short time, even for massive over densities. Accordingly, the 
quintessence effect may be neglected to a first extent. 



26 



From this point on, further effort shall be devoted to the special case, 
w = —1/3, for the following reasons: (i) the overdensity equation of motion 
can analytically be integrated, according to Eq. (69), even in presence of (fully 
clustered) quintessence; (ii) the (matter + quintessence) overdensity behaves 
as a pure matter overdensity within an open CDM universe where r2cDM = 
Qm = 1 ^^(? (e.g., Horellou and Berge, 2005), and the related results may be 
extended to the case under discussion; (iii) overdensities reach turnaround 
and collapse earlier for increasing values of the quintessence equation of state 
parameter, w, and vice versa (e.g., Horellou and Berge, 2005), which makes 
the case under discussion, w = —1/3, a lower limit with regard to the range 
of interest, —1 < w < —1/3, where the maximum difference does not exceed 
about 15% if collapse ends at present for w = —1 (Horellou and Berge, 
2005). Concerning central collapse, the above difference is expected to be 
higly reduced. 

According to Eq. (72), at turnaround a = 1, and the scaling time, tj^^^^, 
coincides with the free-fall time of the matter subsystem, which implies 
Tmax = 1. For the initial configuration, the combination of Eqs. (42), (65a), 
and (70a) yields: 



l + Si 



= 1 - - 
TT 



Aj — Af — arcsin ^1 — A, 



1 - (Qm)-' « Si « l;(74a) 

(74b) 



Using the definition of mean density and matter mass conservation, the 
free-fall time, tmax, via Eq. (65b) takes the expression: 



t 



max — "max 



-1/2 



a. 



-3/2 



(75) 



where (p^)i is the mean density within a sphere of radius, rj, centered on 
a local maximum at recombination epoch. The combination of Eqs. (17), 
(74a), and (75) yields: 



^max(^i) 
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G{phm)i 



(76) 



where {phm)i is the density of the matter Hubble flow at recombination epoch. 
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The mean overdensity, 6, can be related to the initial value, 6i, using mass 
conservation together with Eqs. (19), (70), and (73). The result is: 



1 + 5 / a\ 



l + 5i \ahJ 



3 



- 



(77a) 



ah = ; (77b) 

'"max 

where rh is the radius of a sphere with zero overdensity. 

Starting from the definition of local density within an infinitely thin spher- 
ical shell, and mean density within the volume enclosed by the shell: 

mass conservation within the shell and the volume bounded by the shell, 
with regard to the initial configuration, reads (e.g., Peebles, 1980, Chap. II, 
§ 19C; Ryden and Gunn, 1987; Caimmi, 1990): 



Pm{r) = pmiu) { -) ( — 1 ; (80) 



Pm(r)^Pm{ri){-) ; (81) 



and the combination of Eqs. (19) and (77)-(81) yields: 

where the mean overdensity, S{r), can be determined for a fixed cosmological 
model via Eqs. (69) and (77), and the explicit expression of the partial deriva- 
tive, dr/dri, allows the explicit expression of the local overdensity, 5{r), via 
Eqs. (77). 

The related procedure, involving long but stimulating algebra, has been 
performed in different situations related to CDM cosmologies: overdensity 
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evolution in Einstein-de Sitter universes (Peebles, 1980, Chap. II, § 19C; Ry- 
den and Gunn, 1987; Caimmi, 1989); overdensity evolution in Einstein-de 
Sitter universes with the effects of acquisition of angular momentum on the 
expansion included (Caimmi, 1990); overdensity evolution in high-energy 
(flm ~ 0.1) universes (Andriani and Caimmi, 1991); overdensity evolution in 
high-energy {^1^ ~ 0.1) universes with the effects of acquisition of angular 
momentum on the expansion included (Andriani and Caimmi, 1994). The 
third above mentioned case (Andriani and Caimmi, 1991) shall be considered 
in the current attempt. 

Accordingly, the explicit expression of Eq. (82b) reads: 



where the positive and the negative sign within parentheses are related to 
shell expansion and contraction, respectively. For a formal derivation, see 
Appendix C. 

For density profiles obeying Eq. (5), overdensity central collapse is com- 
pleted at a dimensionless time, T(rpk) = 2, according to Eq. (70b), where rpk 
is the radius of the homogeneous peak at recombination epoch. 

Let Tc(rj) be the dimensionless time related to an infinitely thin spherical 
shell of initial radius, rj, at the end of central collapse i.e. t — 2tmax.{rpk) . 
The combination of Eqs. (65) and (76) produces: 



where (5pk = 5j(rpk) is the overdensity value of the homogeneous peak, and 
6i{ri) is related to the initial density profile via Eqs. (23) and (27). 

The related dimensionless distance, ac{ri), can be determined as the so- 
lution to the transcendental equation deduced from Eq. (70b): 





(84) 



Tc(ri) = 1 - Jacin) - a^(ri) + arcsin 



(85) 
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and the combination of Eqs. (65a) and (74a) yields: 



^c(n) 



Ai(r,) ; 



(86) 



which allows an explicit expression of the radius of an infinitely thin spherical 
shell at the end of central collapse, as: 



rdri) 



ac[ri)ri 
AJrA 



(87) 



where rj is the initial radius at recombination epoch. 

The density profile at the end of central collapse, follows from the partic- 
ularization of Eqs. (82) and (83) to the case of interest. The result is: 



p(r,) l + 5(r,) 1+S,{r,) 
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(88b) 



where global and local overdensities, 5i and 6i, are determined by the knowl- 
edge of the initial density profile at recombination epoch, via Eqs. (23) and 
(24), respectively. 

The overdensity evolution obeying an initial density profile defined by 
Eq. (5), lies between the limiting situations, 6 = and 6 = 3. The special 
case, 6 = 0, corresponds to a Roche sphere i.e. a mass point surrounded 
by a massless atmosphere. Accordingly, rpk = and the hole is present 
from the beginning of evolution. The special case, 6 = 3, corresponds to a 
homogeneous sphere which collapses into a hole after a time, t = 2traa.^{Ri)- 
In the general case, < 6 < 3, the homogeneous peak collapses into a hole 
after a time, t = 2tinax(?'pk), while the underlying envelope virializes due to 
deviations from radial motion. 
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According to the results of Subsect. 2.3, the dimcnsionlcss radius, S = 
R/rp]^, Eq. (8), and the dimensionless densities, Ptr/Ppk and ptr/Ppk, Eqs. (10) 
and (11), respectively, depend only on the input parameters, b and k. For 
changing masses and peak heights, the above mentioned quantities necessar- 
ily remain fixed. With regard to the mass ratio, k, = Mpk/Mtj., an empirical 
relation between hole and hosting spheroid mass (e.g., Ferrarese and Ford, 
2005) cannot be used for for the following reasons. 

First, the above mentioned relation holds for ellipticals and early-type 
spirals, but the disk mass is not taken into consideration. On the other 
hand, the disk is dominant in late-type spirals, and both the proto-bulge 
and the proto-disk, together with the (nonbaryonic) dark matter, must be 
included in the initial density profile at recombination epoch. 

Second, the dark matter is dominant with respect to the baryonic matter, 
and the relation between dark matter and hole mass seems to be nonlinear, 
as: 

^ = 5-810'' (t^^)'^ ; 10 ^ MaWMio ^ 1000 ; (89) 

-'"dark ViUU iViio / / 

where z is the redshift and f3 an exponent to be fixed. Special values are: 
{(3,z) = (0.65,6/29) (Ferrarese, 2002); {(3,z) = (0.27,19/58) (Baes et al., 
2003); and /? = 0.39 (Shankar and Mathur, 2007). In addition, Eq. (89) 
is close to its counterpart cahbrated locally through statistical arguments 
(Shankar et al., 2006) and to its counterparts obtained using semianalytical 
models (Granato et al., 2004; Lapi et al., 2006). It can be seen that a 
redshift-dcpcndent superlinear relation, Mhoie/Miark oc M^^^^, 0.3 ^ f3 ^ 0.7, 
is best suited to represent the luminosity function in active galactic nuclei 
(Shankar et al., 2007). Accordingly, Eq. (89) shall be used in computations. 
An extrapolation to the mass range under consideration, 10~^ < M/Mio < 
10^ via Eq. (89) yields 10"^ < Mhoie/Mio < lO'' in the special case, P = 0.65, 
z = 6/29. 

If the lack of low-mass (lO^'' ~ Mhoic/Mio ~ 10^^) holes is real, a possi- 
ble interpretation could be the following. Secondary perturbations naturally 
arise around a local density maximum (e.g., Ryden and Gunn, 1987) and 
virialize when the hosting overdensity is still expanding. For total clump 
number or typical clump mass independent of overdensity mass, it can be 
seen that the distance between nearest clumps (provided they are uniformly 
distributed within the overdensity volume) increases as the overdensity mass 
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increases, and vice versa. Accordingly, tidal interactions and related de- 
parture from radial motions are expected to be more efficient in low-mass 
(M ~ lO^Mio) overdensities, which could prevent hole formation. 

4 Results 

Computations have been performed up to the end of the central collapse, 
within the mass range, —1 < log(M/Mio) < 6, for mean overdensity heights, 
V = 1,2, 3, 4, and an assumed quintessence equation of state parameter, w = 
— 1/3. The density profile at the end of central collapse, and related values 
of dimensionless times, Tc{ri), and dimensionless radii, adri), are weakly 
dependent on the mean overdensity height, to a major extent as the initial 
shell radius, rj, is selected near the truncation radius. In addition, different 
overdensity heights related to low masses exhibit slightly larger differences 
with respect to their counterparts related to high masses. Also, differences 
in Tc and etc are more pronounced (but still negligible) in comparison to 
differences in density, provided the remaining parameters are kept fixed. 

The ratio of homogeneous peak radius to truncation radius, = rpk/ -Rtr, 
mean density to homogeneous peak density, Ptr/Ppk, and local density at 
truncation radius to homogeneous peak density, ptr/ Ppk, depend only on the 
slope of the external density profile, b, and the ratio of homogeneous peak 
mass to total mass within the truncation radius, k = Mpk/Mtr, according 
to Eqs. (8), (10), and (11), respectively. The dependence within the ranges, 
2.85 < 6 < 3, — 8 < log K. < —2, is shown in Tab. 2. The dependence on the 
power-law density profile exponent, b, for fixed ratio of homogeneous peak 
mass to total mass within the truncation radius, k, is steeper for low k and 
vice versa. 

For proto-galaxies (Af/Mio = 100), slopes b < 2.85 would imply ex- 
ceedingly large initial peak overdensities (5pk > 0.74), and slopes b > 2.996 
would imply exceedingly large cosmic times at the end of central collapse 
{tc > 0.9 Gyr). Mass ratios, outside the above mentioned range would 
be inconsistent with Eq. (89), related to the cases of interest deduced from 
observations and/or supported by theoretical arguments. The homogeneous 
density profile, 6 = 3, makes a useful limiting case. 

In the current model, the hole mass coincides with the homogeneous peak 
mass, and k, can be deduced from Eq. (89). As a reference case, the values 
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Table 2: From up to down of each block: ratio of homogeneous peak radius to 
truncation radius, = rpk/-Rtr (first lines), mean density to homogeneous 
peak density, p^r/Ppk (second lines), local density at truncation radius to 
homogeneous peak density, ptr/Ppk (third hnes), and peak overdensity, 5pk, 
(fourth hnes), for different choices of ratio of homogeneous peak mass to 
total mass within the truncation radius, k = Mpk/M^r, and power- law density 
profile exponent (from top to bottom), b — 2.85 (first block); b — 2.90 (second 
block); b = 2.95 (third block); 6 = 3 (fourth block). 
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Table 3: The homogeneous peak radius, rpk, the ratio of homogeneous peak 
radius to truncation radius, = rpk/-Rtr, and the ratio of homogeneous 
peak mass to total mass within the truncation radius, k — Mpk/Mtr, for 
ovcrdcnsitics with different masses, at recombination epoch, density profile 
defined by Eq. (5), k, deduced from Eq. (89), particularized to the selected 
reference case, b = 2.90 and z) = (0.65, 6/29). The rms overdensity value, 
{SM)i = Sm, and the overdensity radius, Ri = R, are taken from Tab. 1 for 
sake of completeness. 

{(3,z) — (0.65,6/29) shall be chosen, which make Eq. (89) coincide with its 
empirical counterpart deduced from a sample of 37 spiral galaxies (Ferrarese, 
2002), together with the value, b = 2.90 in Eq. (5). The related values of the 
homogeneous peak radius, rpk, and the ratio of the homogeneous peak radius 
to truncation radius, = rpk/i?tr, are listed in Tab. 3 for different masses 
at the beginning of evolution, assumed to be at recombination epoch. Initial 
radii, Ri = R, and rms overdensity values, {SM)i = ^m, are taken from Tab. 1 
for sake of completeness. The above mentioned quantities are independent 
of the overdensity height. 

The homogeneous peak overdensity, 5pk, the local overdensity at the 
truncation radius, Stj-, the homogeneous peak density, Ppk, the overdensity 
turnaround radius, -Rmaxj the cosmic time at turnaround, tmaxi and at the 
end of central collapse, t^, arc listed in Tab. 4 for mean overdensity heights, 
Vi = 1, 2, 3, 4 (from top to bottom in each bloch), and different masses, with 
regard to the selected reference case. The homogeneous peak density, Ppk, 
and the local overdensity at truncation radius, 5tr, decrease only slowly with 
increasing masses, and the same holds for the homogeneous peak overden- 
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Table 4: Homogeneous peak overdensity, (5pk, local overdensity at truncation 
radius, 5tr, homogeneous peak density, Ppk, overdensity turnaround radius, 
-Rmax, cosmic time at turnaround, imax, and at the end of central collapse, 

tc, for mean overdensity heights, Vi = 1,2,3,4 (from top to bottom of each 



sity, 5pk, though to a sUghtly larger extent. In particular, the model density 
profile, expressed by Eq. (5), implies the existence of an underdense external 
region (negative overdensity) up to the truncation radius, at the end of cen- 
tral collapse. Both the overdensity turnaround radius, rmax, and the cosmic 
time at the end of central collapse, t^, decrease for increasing mean overden- 
sity height, z7j. Accordingly, in the light of the model under discussion, hole 
formation occurs within a few hundredths of Gyr, or z > 10. 

The fractional density at the end of central collapse, p/p, as a function 
of the fractional radius, ^ = r / R, at the end of central collapse, for different 
masses and mean overdensity heights, with regard to the selected reference 
case, is plotted in Fig. 1. The density profile of the outer shells, which are 
still expanding together with the universe, exhibits no dependence on the 
overdensity mass. The contrary holds for the inner shells where, for assigned 
fractional radius, the fractional density is an increasing function of the over- 
density mass. The dependence on the mean overdensity height appears to 
be negligible. 

The density, p, as a function of the radius, r, at the end of central collapse, 
for different masses and mean overdensity heights, with regard to the selected 
reference case, is plotted on a logarithmic plane in Fig. 2. The density profile 
is universal, with the exception of the inner shells, where a higher density 
is attained. The dependence on the mean overdensity height appears to be 
negligible. 

The dimensionless time, r = Tc(ri) = tc/tma^{ri), as a function of the 
fractional radius, = r / R, at the end of central collapse, for different masses 
and mean overdensity heights, with regard to the selected reference case, is 
plotted in Fig. 3. The dimensionless time related to the outer shells, which are 
still expanding together with the universe (r ~ 0), exhibits no dependence on 
overdensity mass. The contrary holds for the inner shells (r ~ 0.5) where, for 
assigned fractional radius, the dimensionless time is an increasing function 
of the overdensity mass. In other words, with respect to a selected fractional 
radius, the homogeneous peak collapses "faster" in low-mass overdcntity than 
in large-mass overdensities. The dependence on the mean overdensity height 
appears to be negligible. 

The dimensionless distance, a — CKc{ri) — r{tc)/rmaxiri), as a function of 
the fractional radius, ^ = r / R, at the end of central collapse, for different 
masses and mean overdensity heights, with regard to the selected reference 
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Figure 1: The fractional density, p{r)/p — p{r)/ < p >, as a function of the 
fractional radius, ^ = r/i?, at the end of central collapse, for different masses 
and mean overdensity heights, z7j = u, with regard to the selected reference 
case. Different symbols correspond to the following values of log(M/Mio): 
— 1 (squares); (triangles); 1 (asterisks); 2 (Greek crosses); 3 (diamonds); 4 
(St. Andrew's crosses); 5 (reversed triangles); 6 (dots). 
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Figure 2: The decimal logarithm of the density, log[p/(Mio/kpc^)], as a 
function of the decimal logarithm of the radius, log(r/kpc), at the end of 
central collapse, for different masses and mean overdensity heights, Vi = u, 
with regard to the selected reference case. Caption of symbols as in Fig. 1. 
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Figure 3: The dimensionless time, r = rc(rj) = tc/^max( 

function of 

the fractional radius, C, = r/R, at the end of central collapse, for different 
masses and mean overdensity heights, with regard to the selected reference 
case. Caption of symbols as in Fig. 1. 
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Figure 4: The dimensionless distance, a = a^ri) = r(tc)Mnax( 
function of the fractional radius, ^ = r/R, at the end of central collapse, for 
different masses and mean overdensity heights, with regard to the selected 
reference case. Caption of symbols as in Fig. 1. 



case, is plotted in Fig. 4. The dimensionless distance related to the outer 
shells, which are still expanding together with the universe (q;~ 0), is a de- 
creasing function of the overdensity mass, according to the initial conditions 
via Eq. (86). The above mentioned effect is less evident moving inwards, un- 
til no appreciable mass dependence is shown. The trend is reversed for the 
still expanding inner shells where the dimensionless distance is an increasing 
function of the overdensity mass. For fixed fractional radius, r{tc)/ R{tc), 
the inner shells are closer to maximum expansion, or end of contraction, 
for large-mass overdensities with respect to low-mass overdensities. In other 
words, for shells with fixed fractional radius, the homogeneous peak collapses 
"faster" in low-mass overdensities with respect to large-mass overdensities: 
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at the end of central collapse, shells which are turning around correspond to 
a smaller initial fractional radius for low-mass overdensities, with respect to 
large-mass overdensities. The fractional radius at the end of central collapse, 
^ = r/R, where the fractional distance, a, shows no appreciable dependence 
on the overdensity mass, is a decreasing function of the mean overdensity 
height, passing from ^ ~ 0.9 for 77^ = 1 to ~ 0.7 for I7j = 4. 

The effect of a change of the model parameters on the results, is shown 
in Fig. 5 for the dimensionless time, r = Tc{ri) — ic/*max(n); ^ a function of 
the fractional radius, ^ = r/R, at the end of central collapse, and in Fig. 6 
for the dimensionless distance, a = ac(^j) = i^{tc)/i^niax{i^i), as a function of 
the fractional radius, ^ = r/R, at the end of central collapse, respectively, 
for different masses. The selected reference case is placed on the middle 
right, which is the same as in Fig. 3 and Fig. 4 (bottom right), respectively. 
The changes in parameter values with respect to the selected reference case, 
are listed on each panel. 

The limiting situation of homogeneous initial density profiles, 6 = 3, is 
shown for comparison only, as it would imply a cosmic time at the end of 
central collapse, longer than the age of the universe for large-mass overden- 
sities (M/Mio > 10^). Slopes within a fiducial range, 2.85 < b < 2.95, make 
little change on both the dimensionless time and the dimensionless distance, 
as shown in Figs. 5 and 6. A weak trend is exhibited, where the evolution 
of different shells is closer in shallower density profiles and vice versa, as 
expected. 

According to the assumed hole-hosting dark matter halo mass relation, 
Eq. (89), the hole mass is larger (for assigned hosting dark matter halo mass) 
for larger exponent, f3, and/or rcdshift, z, and vice versa, which implies closer 
evolution of different shells for low-mass holes and vice versa. An inspection 
of Figs. 5 and 6 shows that it is the case, and the related trend is more visible 
than in the case of changing slope of the initial density profile, b. 

In the light of the model presented in the current attempt, holes are the 
first structures which form in proto-galaxies, when the upper shells are still 
collapsing or expanding. 
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Figure 5: The dimensionless time, r = Tc(rj) = tc/^max(?^i), as a function of 
the fractional radius, ^ = r / R, at the end of central collapse, for different 
masses, mean overdensity height v — V — 3, and different slopes of the initial 
density profile, b, Eq. (5), exponent, /3, and redshift, z, appearing in the 
assumed hole-hosting dark matter halo mass relation, Eq. (89). The selected 
reference case is placed on the middle right (same as in Fig. 3, bottom right). 
The parameter values which have been changed with respect to the reference 
case, are listed on each panel. The limiting situation of homogeneous initial 
density profiles, 6 = 3, is shown for comparison only. Caption of symbols as 
in Fig. 1. 
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Figure 6: The dimensionless distance, a = adri) = r(tc)/Tmax(?^i), as a func- 
tion of the fractional radius, ^ = r/R, at the end of central collapse, for 
different masses, mean overdensity height u — V — 3, and different slopes of 
the initial density profile, b, Eq. (5), exponent, /3, and redshift, z, appearing 
in the assumed hole-hosting dark matter halo mass relation, Eq. (89). The 
selected reference case is placed on the middle right (same as in Fig. 4, bot- 
tom right). The parameter values which have been changed with respect to 
the reference case, are listed on each panel. The limiting situation of homo- 
geneous initial density profiles, 6 = 3, is shown for comparison only. Caption 
of symbols as in Fig. 1. 
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5 Discussion 



Gravitational collapse may be a viable mechanism for hole generation, due 
to two orders of reasons. First, remnants of sufficiently massive Type II su- 
pernovae are thought to be black holes (e.g., Nomoto et al. 2005; Ohkube 
et al. 2006). Second, the mean density of a (spherical-symmetric) black hole 
is a decreasing function of the gravitational radius. Low-mass (m ~ 2Om0) 
Type II supernovae leave a neutron (or more exotic) star remnant, which 
is sustained by the Fermi pressure within a radius about three times larger 
than the gravitational radius. The Fermi pressure makes the collapsing core 
bounce at a density about twice the nuclear density, and later attain a viri- 
alized configuration. On the other hand, the bounce would occur inside the 
gravitational radius for more massive (m ~ 20-25mo) Type II supernovae, 
which implies black hole formation. If (nonbaryonic) dark matter is mainly 
made of fcrmions, a similar mechanism could hold. Low-mass (M ~ lOMio) 
dark matter overdensities could leave a fermion ball remnant at the end of 
central collapse, while larger masses could imply hole generation (e.g., Viol- 
lier 1994; Munyaneza and VioUier 2002; Munyaneza and Biermann 2006). 

Radial motions are better preserved in coUisional fluids where energy dis- 
sipation takes place, such as in supernovae, than in collisionlcss fluids where 
no energy dissipation occurs, such as dark matter overdensities. Clump for- 
mation and tidal interactions between clumps or from neighbourhing over- 
densities, convert ordered radial motions into random orbital motions. On 
the other hand, the local maximum at the overdensity deep interior, is a 
very special place where both tidal interactions and acquisition of angular 
momentum may safely be neglected. Accordingly, radial motions could be 
preserved up to hole formation. 

The extrapolation of the hole-hosting dark matter halo mass relation, 
Eq. (89), to early cosmic times, has been considered as the most viable ex- 
pression of the homogeneous peak to overdensity mass ratio. The related 
value could be overstimated if only the inner part of the homogeneous peak 
is the hole progenitor, and understimated if substantial mass accretion onto 
the hole works. The extrapolation of Eq. (89) in the above mentioned sense, 
implies the mass of the hole progenitor is proportional to the overdensity 
mass, at flxed redshift. If the total clump number or the typical clump mass 
does not appreciably depend on the overdensity mass, clumps within low- 
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mass ovcrdensities are expected to be efficient in exerting tidal actions and 
related departure from radial motions. Accordingly, central collapse and hole 
formation could be inhibited, in alternative to the formation of a fermion ball. 

Though the assumed density profile at the beginning of evolution, Eq. (5), 
has no physical motivation, still it cannot be considered as completely ad hoc. 
In fact, it is lying between two limiting situations, related to a homogeneous 
and a generalized Roche sphere (a homogeneous sphere surrounded by a 
massless atmosphere), respectively. It is defined by two parameters: a power- 
law exponent, b, and a fractional mass, k, = Mpk/Mtr. The special cases, 
b = and 6 = 3, are related to generalized Roche and homogeneous spheres, 
respectively. A fiducial range, 2.85 < b < 2.95, has been chosen for the 
following reasons: low values (6 < 2.85) would imply exceedingly high peak 
overdensities, while large values {b > 2.95) would imply exceedeingly long 
cosmic times at the end of central collapse. The main feature is the presence 
of a homogeneous peak which turns around and collapses into a black hole, 
and an inhomogeneous envelope which, after expansion, collapses and relaxes. 

The fractional mass, k, has been selected as a function of both the over- 
density mass and the redshift, expressed by Eq. (89), provided overdensity 
and homogeneous peak mass coincide with dark halo and hole mass, re- 
spectively. Recent empirical correlations are related to special choices of 
the parameters in Eq. (89), namely: {f3,z) = (0.65,6/29) (Fcrrarcsc, 2002); 
{(3,z) = (0.27,19/58) (Baes et al., 2003); (3 = 0.39 (Shankar and Mathar, 
2007). Also, Eq. (89) is close to its counterpart calibrated locally through 
statistical arguments (Shankar et al., 2006) and to its counterpart obtained 
using semianalytical models (Granato et al., 2004; Lapi et al., 2006). In con- 
clusion, the assumed density profile at the beginning of evolution, Eq. (5), 
may be considered phenomenological one. 

Strictly speaking, the model discussed in the current attempt holds for 
isolated overdensities in a secondary infall scenario, where the inner and 
denser regions first virialize while the outer and less dense regions are still 
expanding (Gunn 1977). On the other hand, the formation of dark haloes 
is characterized by accretion of smaller subunits or merger between systems 
of comparable size, and identification of dark matter haloes in simulations 
appears problematic (e.g., Bett et al. 2007). Unexpectedly, a recent inves- 
tigation shows that the secondary infall scenario provides a valid theoretical 
framework for calculating the structure and evolution of dark matter haloes 
in an expanding universe, and that its predictions with respect to the density 
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profiles are in close agreement with full iV-body simulations (Ascasibar et al. 
2007). In this view, most of the diversity in density profiles is contributed 
by the scatter in their primordial counterparts rather the scatter in angular 
momentum (Ascasibar et al. 2007) i.e. the nature of the initial conditions 
rather than the details of the evolution. 

Though hole gravitational radius cannot be resolved by numerical simu- 
lations, the resulting dark halo density profile appears to be cusped at the 
centre (e.g., Navarro et al. 1995, 1996; Moore et al. 1998, 1999; Fukushige 
and Makino 2001, 2003; Diemand et al. 2004; Reed et al. 2005), but recent 
results also allow cored density profiles (e.g., Fukushige and Makino 2004; 
Navarro et al. 2004; Merritt et al. 2005). If the presence of a central cusp 
in simulated dark matter haloes is a real effect instead of an artefact due to 
computer codes, a physical interpretation could be hole formation at the end 
of central collapse. 

Central collapse within local density maxima is one way of avoiding the 
problem of haloes at high redshifts {z ~ 6) that have had insufficient time 
to grow at the Salpeter rate from solar or intermediate mass black hole pro- 
genitors (e.g., Willott et al., 2003; Miller et al., 2006). In this view, holes 
are the first structures formed in evolving overdensities, while proto-haloes, 
proto-bulges, and proto-disks are still relaxing. The active (quasar) phase 
begins just after hole formation, when the outstanding shells are still collaps- 
ing, and may occasionally be enhanced by merger events. The cosmic time 
at the end of central collapse does not exceed a few hundredths of Gyr, see 
Tab. 4, in agreement with observations of high-redshift quasars. Virialization 
takes place later, at a cosmic time equal to about 3 Gyr, or z ^ 2, provided 
M/Mio = 500 [deduced from Cromm et al., 2005, using h = 0.65 in accor- 
dance with Eq. (35)] and u = 3 are typical for luminous quasars, see again 
Tab. 4. The occurrence of virialization, gas exhaustion due to star formation, 
and absence of merger events in the hosting spheroid component, mark the 
beginning of the quiescent phase. 

The simple model used in the current attempt deals with isolated overden- 
sities where neither accretion nor merging occurs. Hole formation via dark 
matter collapse implies a limited range of mass, where the lower limit is re- 
lated to the stability of a fermion ball (provided dark matter is mainly made 
of fermions) or to the occurrence of tidal effects between clumps, and the 
upper limit is related to the end of central collapse at a cosmic time equal 
to the age of the universe. Accordingly, very massive holes are expected 
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within clusters of galaxies and superclusters. Though the model is limited to 
hole formation, still some qualitative considerations on the evolution can be 
performed. 

The picture of holes forming at high redshifts and remaining largely un- 
changed since then, is not consistent with the observed luminosity density 
of active galactic nuclei and the inferred local hole density, unless the radia- 
tive efficiency of the luminosity of active galactic nuclei is unfeasibly high 
(Marconi et al., 2004). So it appears that holes have been continuing to ac- 
crete mass during the cosmic epochs in which dark matter haloes and their 
associated galaxies have also continued to be assembled. 

In the light of the model used in the current paper, the above scenario 
still holds with regard to the baryonic precursor of the hole mass budget. 
On the other hand, little change occurs since hole formation with regard 
to the (dominant) nonbaryonic precursor. In fact, the basic assumption is 
that seed holes arc created at the end of central collapse (z > 10) instead 
of in major mergers which, in addition, trigger an episode of gas accretion 
in overdensities with a pre-existing hole. In other words, mass accretion 
occurs after hole formation, while in earlier attempts mass accretion occurs 
during hole formation (e.g., Di Matteo et al., 2003; Bromley et al., 2004; 
Miller et al., 2006; Lapi et al., 2006; Li et al., 2007). Baryon accretion may 
highly be reduced by central starburst or radiative feedback (e.g., Ciotti and 
Ostriker, 2007). Nonbaryonic matter accretion may not significantly increase 
hole mass, due to the lack of a mechanism to dissipate angular momentum, 
and it may be important only in early times, 2; ~ 30 (Mack et al., 2007). 

Several previous models have considered the possibility that holes form 
from low-mass black holes related to pop. II (e.g., Haiman and Loeb, 2001) 
or pop. Ill (e.g., Volonteri et al., 2003) stars. However, there are vari- 
ous problems associated with the above mentioned scenarios (e.g., Haehnelt, 
2003; Pelupessy et al., 2007). More specifically, some mechanism is required 
to facilitate the migration of seed black holes to the centre of their hosting 
spheroid. 

Another class of models focused on the possibility that holes form directly 
from the collapse of a large gaseous cloud placed in the central overdensity 
region (e.g., Rees, 1984; Haehnalt and Rees, 1993; Silk and Rees, 1998; 
Bromley et al., 2004), where the problem is the need to avoid fragmentation 
and related departure from radial motions during the collapse. In particular, 
it is assumed (Bromley et al., 2004) that the collapse is also connected to 
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overdensity major mergers, in the sense that the same tidally stripped gas 
which falls into the centre of a post-merger overdensity and can fuel quasar 
accretion, is also viewed as the source of the initial hole formation. 

The model used in the current paper follows the same line of thought 
but is different to many respects, namely: (i) hole formation is (mainly) 
due to nonbaryonic matter instead of gas; (ii) the initial driving mechanism 
is contraction of local density maxima after turnaround instead of major 
overdensity merger; (iii) quasar accretion is initially fuelled by gas secondary 
infall instead of tidally stripped gas during overdensity major mergers. 

As reported in an earlier attempt (Bromley et al., 2004), any potential 
model of quasar and active galactic nucleus formation finds itself faced with 
three major unknowns. First: where, how and with what mass do the initial 
holes form? Secondly: what events trigger their subsequent fuelling and 
growth, how much fuel do they supply and how efficiently is it converted into 
radiative energy? Thirdly: how does feedback (from star formation or the 
active galactic nucleus activity itself) regulate and possibly even check their 
growth? 

The current paper answers the first point raised above: hole formation 
is the result of central collapse in local density maxima, mainly due to non- 
baryonic dark matter. Subsequent fueling, growth, and feedback, are due 
to baryonic matter, the sole which can radiate and make quasars and active 
galactic nuclei shining. On the contrary, the hole mass budget is mainly 
determined by its nonbaryonic precursor. In this view, current models on 
active galactic nuclei formation and co-evolution with hosting galaxies, are 
expected still to hold with regard to the second and the third point raised 
above. 

6 Conclusion 

Unsustained matter distributions unescapely collapse unless fragmentation 
and centrifugal or pressure support take place. Starting from the above ev- 
idence, supermassive compact objects at the centre of large-mass galaxies 
(defined as "holes") have been conceived as the end-product of the gravita- 
tional collapse of local density maxima (defined as "central collapse") around 
which positive density perturbations (overdensities) are located. At the be- 
ginning of evolution, assumed to occur at recombination epoch, local density 
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maxima have been idealized as homogeneous peaks, while the surrounding 
envelopes have been described by a power-law density profile, p(r) oc r'^~^, 
< 6 < 3, where 6 = represents a massless atmosphere and 6 = 3 a 
homogeneous layer. The dependence of the density profile on a second pa- 
rameter, chosen to be the ratio between peak and total (truncated) mass, 
K = Afpk/Aftr, has been analysed. 

Overdensity evolution has been discussed in the context of quintessence 
cosmological models, which should be useful in dealing with the virialized 
phase. Aiming to describe the central collapse, further investigation has 
been devoted to a special case where the quintessence effect is equivalent to 
additional curvature {w = —1/3), and overdensities exhibit the selected den- 
sity profile at recombination epoch. A redshift-dependent, power-law relation 
between hole and (nonbaryonic) dark halo mass has been used to express the 
dependence of the fractional mass, k = Mpk/Mtr, on the overdensity mass, 
M = Mtr, where the homogeneous peak and overdensity mass are related to 
the hole and dark halo mass, respectively. 

Computations have been performed for a wide range of masses, — 1 < 
log(M/Mio) < 6, and mean overdensity heights, 1 < i^j < 4, up to the end of 
central collapse, for the following quantities: homogeneous peak overdensity, 
5pk; turnaround radius, rmax! cosmic time at the end of central collapse, 
tc, density profile, p[r{tc)]; logarithmic density profile, log[p/(Mio/kpc'^)]; 
ratio of cosmic time at the end of central collapse to turnaround time of the 
related shell, Tc = tc/^l^maxl^'i)]; ratio of shell radius at the end of central 
collapse to shell radius at turnaround, etc = r{tc) / rmax{ri) ■ With regard to 
the last two above mentioned quantities, additional computations have been 
performed for different slopes of the initial density profile, b, and different 
parameters appearing in the hole-hosting dark matter halo mass relation, (3 
and z, Eq. (89). 

The central collapse has been found to end in early times, no longer than 

a few hundredths of Gyr, which implies hole formation when proto-haloes, 
proto-bulges, and proto-disks are still relaxing. No appreciable change has 
been found in the evolution (up to the end of central collapse) of differ- 
ent mean overdensity heights related to equal masses. On the other hand, 
it has been recognized that homogeneous peaks collapse (in dimensionless 
coordinates) "faster" with respect to surroundings envelopes, in low-mass 
overdensities than in large-mass overdensities. In conclusion, it has been in- 
ferred that gravitational collapse of homogeneous peaks within overdensities 
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may be a viable mechanism for hole generation. 
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Appendix 



A On the light velocity in vacuum 

Strictly speaking, light travels across the vacuum (intended in philosophi- 
cal sense, as absence of everything) only in classical (Newton) and special- 
relativistic (Einstein) mechanics, where it is conceived as made of corpuscles. 
If, on the other hand, hght is conceived as waves (Huygens), then it propa- 
gates across a homogeneous medium called ether. In general-relativistic me- 
chanics (Einstein) ether is replaced by gravitational field of the universe as a 
whole. In quantum mcchanichs (Bohr, de Broglie, Heisemberg, Schrodingcr) 
light can be considered as made of either corpuscles or waves, and the prop- 
agation takes place across the quantum void, which hosts some form of en- 
ergy, and then has not to be intended in the above specified philosophical 
sense. Light velocity is maximum in vacuum (or ether, or quantum void), 
while it decreases for increasing refraction index when propagation occurs 
within a (homogeneous) medium, due to photon interaction with medium 
constituents. 

Light velocity has been deduced from experiments on the Earth or within 
the solar system at most. According to current QCDM cosmologies, the solar 
system is embedded within (nonbaryonic) dark matter and dark energy. The 
current value of the light velocity coincides with the maximum attainable 
value only if no interaction is assumed to occur between (i) photons and 
dark matter, and (ii) photons and dark energy. If otherwise, light velocity 
could be higher in intergalactic voids and/or in early (cosmic) times, where 
the effect of the dark energy was negligible. Even in pure matter universes, 
light propagation takes place across the quantum void, where some form of 
energy is present. 

In conclusion, the term "light velocity in vacuum" has to be intended 
as "light velocity in absence of baryonic matter" under the assumption that 
photons interact with neither (nonbaryonic) dark matter, nor dark energy. 

B The [kpc Mio Gyr] system of measure 

In dealing with large-scale celestial objects, such as galaxies or cluster of 
galaxies, it may be convenient to use an appropriate system of measure in- 
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stead of the standard astrophysical one, [cm g s]. To this aim, unit length, 
unit mass, and unit time, shall be taken as: 

1 kpc = 3.085 677 580 666 31 10^^ cm ; (90) 
Mio = 10^°mQ = 1.989 1 10^^ g ; (91) 
lGyr = 10^ y = 3.153 6 lO^S ; (92) 

and the related "mathematical" year is defined as 1 y = (365 x 24 x 60 x 60)s, 
independent of the Earth revolution around the Sun. 
According to Eq. (90), the astronomical unit is: 

1 AU = 1.495 978 706 91 10^^ cm = 4.848 136 812 10 lO"'^ kpc ; (93) 

using Eqs. (90) and (92), the velocity unit reads: 

1 kpc Gyr"^ = 0.978 461942118 946 6 km s"^ ; (94a) 
1 km s"^ = 1.022 012 156 992 443 kpc Gyr"^ ; (94b) 

and the light velocity in vacuum is: 

c = 2.997 924 58 10^ km s~^ = 3.063 915 37 10^ kpc Gyr"^ ; (95) 

on the other hand, the value of the gravitation constant is deduced from 
Eqs. (90)-(92), as: 

G = (6.672 59 ^ 0.000 85) g"^ cm^ s'^ = (44927.5 ^ 5.7) M^^ kpc^ Gyr"^ . 

(96) 

Using Eqs. (90) and (91), the unit density reads: 

1 Mio kpc"^ = 6.770 3 lO'^^ g cm"^ ; (97a) 
lgcm-3 = 1.4770 10^^ Miokpc"^ ; (97b) 

and the value of the mean solar density is: 

Pq = 1.411 gcm-=^ = 2.083 10^^ Mio kpc~^ ; (98) 

where Rq = 6.96 lO^'' cm has been assumed. 

Using Eqs. (90)-(94), the unit energy and the unit specific energy read: 

1 Mio kpc^ Gyr-^ = 1 .904 3 10^° g cm^ s"^ ; (99a) 
1 g cm^ s-2 = 5.251 2 10"" Mio kpc^ Gyr"^ ; (99b) 
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= 9.573 877 719 424114 10^ cm^s"^ ; (100a) 
1 cm^ = 1.044 508 849 294 298 10"^° kpc^ Gyr-^ ; (100b) 

similarly, the unit angular momentum and specific angular momentum read: 

1 Mio kpc^ Gyr^^ = 6.005 5 10^^ g cm^ s"^ ; (101a) 
1 g cm^ s"^ = 1.665 1 10~™ Mio kpc^ Gyr"^ ; (101b) 

1 kpc^ Gyr"^ = 3.019 218 077 027 732 10^^ cm^ ; (102a) 
Icm^s"^ = 3.312115 834 257490 10"^^ kpc^Gyr"^ ; (102b) 

accordingly, velocities, energies, and angular momenta, may be translated 
from [cmgs] to [kpcMioGyr] system of measure and vice versa. 

C Evolution of bound overdensities in CDM 
universes 

Overdensity evolution in CDM universes is determined by Eqs. (17), (61), 
(65), (70), and (77)- (82). To perform calculations, an explicit expression of 
the function, g{r), defined by Eq. (82b), is needed. To help the reader, the 
procedure used in earher attempts (Andriani and Caimmi, 1991, 1994) shall 
be repeated here for the case of interest i.e. bound perturbations in an open 
universe. 

To this aim, it is useful to define the dimensionless variables (Andriani 
and Caimmi, 1991): 

; ^'=7tV ; (103a) 

n {tffh 

(the index, h, denotes zero density excess) which are related to their coun- 
terparts, expressed by Eqs. (65), via Eq. (74) as: 

a^A,a' ; r = (1 + ^O'^'Ap^' ; (104) 
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with regard to an infinitely thin spherical shell of initial radius, rj. 

In addition, both the cosmic time, t, and the free-fall time at the beginning 
of evolution, {tff)i, via Eqs. (74) are independent of rj, which implies null first 
derivatives: 

dr' drl 

dr-i dr-i ^ ^ 
Using Eqs. (103), (104), and the trigonometric identity: 

1 



- arcsin(l — 2x) + — = arcsin Vl — x 
2 4 

the following expression of Eq. (69) holds: 



< X < 1 ; 



T - T, 



T-(l + -^A,a'-(A,a')2+ ^A^ - Af 

TT L 

-|-^ arccos(l — 2AjQ;') — ^ arccos(l — 2Aj) 

where the plus is related to expansion, and the minus to contraction. 
The combination of Eqs. (82b) and (103a) yields: 



(106) 



gyr) 



a 



d{a'ri) 



1 + 



a' dn 



in terms of the dimensionless radius, a' . 

For spherical-symmetric overdensities, Eq. (20) reads: 

5 — ^ [ S(r)r'^ dr ; 
Jo 

and a derivation with respect to r produces: 



(107) 



(108) 



dr 



(109) 



Keeping in mind that the density parameter, {flm)i, is independent of r^, 
a derivation with respect to on both sides of Eq. (42b) yields: 



5A, dai 



dri 



dri 



-1 + A, d5i 
l+'Si dri 



(110) 



59 



owing to Eqs. (42) and (65a), where the former has been generahzed to the 
sphere bounded by a generic isopycnic surface of initial radius, r^. 

Using Eqs. (103)-(105) and (107)-(110), a derivation with respect to 
on both sides of Eq. (106), after a lot of algebra, produces the following 
expression of g{r): 



r . 1 1-A, 5, ( 5,\\ Aj TjAia'-iAia')^ 
[gir)]-' = 1 + 3 , \ ' 1 - =i 1 - ^ ^ 



X 



(111a) 



A' \ J V 3 1 - Ai 

[giriT'^i ; (111b) 

for further details and exhaustive investigation on open, fiat, and closed CDM 
universes hosting bound, zero-energy, and unbound overdensities, refer to 
Andriani and Caimmi (1991). 

Turning again to the earlier dimensionless coordinates, a and r, the com- 
bination of Eqs. (65), (104), and (111) yields Eq. (83). 
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